Growth of Bose-Einstein Condensates from Thermal Vapor 
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We report on a quantitative study of the growth process of 
*''Rb Bose-Einstein condensates. By continuous evaporative 
coohng we directly control the thermal cloud from which the 
condensate grows. We compare the experimental data with 
the results of a theoretical model based on quantum kinetic 
theory. We find quantitative agreement with theory for the 
situation of strong cooling, whereas in the weak cooling regime 
a distinctly different behaviour is found in the experiment. 
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The non-equilibrium path to the formation of a Bose- 
Einstein condensate is one of the most intriguing phys- 
ical processes in ultra cold quantum gases. When the 
critical temperature of this phase transition is crossed, 
the atoms accumulate in the ground state of the system, 
which leads to the characteristic dense core in the atomic 
distribution and to the long-range phase coherence. It 
is an experimental challenge to monitor this process in 
time, aiming for a precise understanding of the forma- 
tion of Bose-Einstein condensates. Experimental studies 
of the growth of Bose-Einstein condensates have so far 
either focused on verifying the effect of bosonic stimula- 
tion 1^ or on the case of attractive interactions between 
the atoms, which gives rise to special dynamics [|| . How- 
ever, the key role of the thermal cloud, from which the 
condensate grows, has not yet been addressed in experi- 
ments. 

In a pioneering theoretical study Kagan et al. Q pro- 
posed the appearance of a quasi-condensate with strong 
phase fluctuations prior to the formation of long-range 
phase coherence. Stoof |^ presented the idea of nucle- 
ation followed by a growth according to a kinetic equa- 
tion. A simple quantitative theory has been developed 
by Gardiner et al. on the basis of quantum kinetic 
theory. The key equation describing the time dependence 
of the number A'o of atoms in the condensate reads 



No = 2W- 



-{No) { (l - e[^c(^«)-^l/'^«^) A^o + l} , (1) 



where /ic(Ao) is the chemical potential of the conden- 
sate, /i the chemical potential of the thermal vapor, and 
W'^{No) an analytically known rate factor [Q. Starting 
from TVo = the condensate starts to grow at a finite 



rate, which turns quickly into exponential growth, before 
it reaches its asymptotic equilibrium value when the con- 
densate chemical potential approaches that of the ther- 
mal cloud. This equation describes only the growth of the 
population in the ground state of the system and it ne- 
glects certain scattering processes among atoms in higher 
lying states. Especially in the initial period the redistri- 
bution of population between higher lying states becomes 
important and modifies the growth Q. When the scat- 
tering terms are included one finds reasonable agreement 
with the experiment described in Hi, but qualitative and 
quantitative discrepancies remain 0,^. Moreover, equa- 
tion neglects the dynamics of the thermal component 
by assuming its chemical potential to be stationary. The 
infiuence of the dynamics of the thermal cloud was in- 
cluded in the recent models |^J^. A purely numerical 
approach by Monte Carlo simulations is presented in . 

We have identified and studied different stages of for- 
mation when the Bose-Einstein condensate grows from 
a thermal gas. In the experiment a cloud of thermal 
atoms is prepared above the critical temperature and 
then cooled into the quantum degenerate regime by con- 
tinuous evaporation for a duration of up to 6 s. This 
strategy allows us to observe the formation process in 
slow motion. A distinct growth characteristic is revealed 
by directly controlling the properties of the thermal va- 
por. In comparison, previous experimental and most the- 
oretical work addressed the situation in which a short 
radio-frequency pulse is applied to remove hot atoms, 
then turned off leaving the system to evolve toward ther- 
mal equilibrium without any further forced evaporation. 

For a quantitative understanding it is crucial to pre- 
pare a well-defined initial state of the thermal cloud. As 
described previously |l^, we load about 10^ atoms in the 
\F — \,mp = —1) hyperfine ground state into a mag- 
netic Quadrupole-Ioffc configuration trap (QUIC-trap). 
The trapping frequencies are luj_ = 2tt x 110 Hz in the 
radial and uj^ — 2n x 14 Hz in the axial direction. In the 
magnetic trap we perform forced evaporative cooling for 
about 25 s by continuously lowering an RF frequency to 
a final value of lynpo— 2120kHz. This prepares a purely 
thermal cloud of A'i=(4.2±0.2) x 10^ atoms at a temper- 
ature of ri=(640±30) nK. The minimum of the trap is 
determined by atom laser output coupling ||l3] and corre- 
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spends to a radio-frequency of 1955 kHz. The stability of 
the magnetic trap enables us to approach the formation 
of the condensate very slowly. The measured stability 
psf and the good reproducibility of the experiment show 
that the formation of the condensate is not triggered by 
fluctuations of the magnetic field. When operating the 
trap at full current, we have measured a slow upwards 
drift of the bottom of the trapping potential of 5kHz/s, 
which is due to thermally induced motion of the trap 
coils. 

After this preparation of a purely thermal cloud, the 
atoms are cooled through the BEC phase transition. The 
radio-frequency field jumps abruptly to a lower frequency 
J^RF,! and remains at this value for a period of up to 6 
seconds. This limits the depth of the magnetic trap to 
Ccut) which is expressed in units of Ti by the evaporation 
parameter rj = ecut/ksTi. Atoms with an energy larger 
than Cent are lost from the trap due to spin flips into a 
non-trapped state. Equilibration of the cloud by elas- 
tic collisions leads to a decrease of temperature and the 
Bose-Einstein condensate forms during this process. 

The time evolution of the condensation is probed by 
releasing the atoms from the magnetic trap after different 
evolution times, for otherwise identical repetitions of the 
experiment. The atom cloud is imaged by absorption 
imaging after 17 ms of ballistic expansion. Atom num- 
bers are extracted from the absorption images by two- 
dimensional fits to the atomic density distribution. The 
thermal cloud is fitted by a Bose-Einstein distribution 
under the assumption of zero chemical potential, and the 
condensed fraction is fitted by a Thomas- Fermi distribu- 
tion. The error in the determination of the atom number 
is estimated to be below 10%. Here, the "condensate 
number" refers to the integral under the Thomas-Fermi 
profile. Temperatures are determined from a Gaussian 
fit to the far out wings of the thermal atom cloud. The 
statistical error on the temperature determination is 5%, 
whereas the systematic error is 10-15%. 

The initially prepared atomic sample corresponds to 
77 = 6, where the gravitational sag of the atom cloud must 
be taken into account. Gravity pulls the atomic cloud 
away from the center of the magnetic trap and thereby 
reduces the dimensionality and the efficiency of evapo- 
rative cooling. We have performed experiments in the 
range 0.75 < < 4.8, corresponding to RF frequencies 
between 1990 kHz and 2095 kHz. For i] > 4.6 we could 
not obtain Bose-Einstein-condensates within the 6 s pe- 
riod of the measurement due to the low efficiency of the 
cooling. For rj — 0.75 we did not observe the formation 
of a condensate since the atom number of the initial ther- 
mal cloud was reduced too much for evaporative cooling 
to proceed. 

Figure ^ shows the measured data for an evaporation 
parameter rj = 1.4. Every data point is an average 
over three measurements and the error bars show the 
statistical error. The final temperature of the cloud is 



(220±20) nK and the final condensate fraction is (55 ± 
7)%. The solid line is a numerical calculation of the con- 
densate growth according to a model based on quantum 
kinetic theory, where the infiuence of the thermal atoms 
is taken into account The model was improved to 
include the continuous evaporation of atoms, the gravi- 
tational sag of the atomic cloud, the experimentally ob- 
served magnetic field drift, and three-body loss of atoms 
from the condensate jl^. We find quantitative agree- 
ment between the model and the measured data with- 
out any free parameters entering into the calculation, 
which suggests that the model describes this regime of 
the formation of Bose-Einstein condensates sufficiently. 
The best agreement between data and theory has been 
obtained for the calculation with the starting conditions 
Ni = 4.4 X 10^ and T, = 610 nK or Ni = 4.2 x 10^ and 
Ti = 610 nK, both of which lie within the range of the 
experimental uncertainties. 

We have compared the time of the onset of Bose- 
Einstein condensation after the beginning of the final 
radio frequency cooling with the calculated curves and 
find excellent agreement for the situation of strong cool- 
ing (see Fig. H). To quantify the initiation time ti for the 
measured data we have fitted the function 



Noit) = Noj ■ (l - e-r(t-*.)) 



(2) 



which describes the final stage of the condensate growth 
towards equilibrium (see Fig. F describes a relax- 
ation rate towards equilibrium. This functional form was 
suggested in Ref. ||l^ and subsequently found to be a 
reasonable approximation to more sophisticated calcula- 
tions 10,^. To account for the range of validity of the 
equation we have restricted the fitting range to conden- 
sate numbers larger than 30% of Nqj, where we have 
obtained good agreement with the data. For low values 
of the evaporation parameter, i.e. rapid cooling, the ini- 
tiation time approaches a minimum of 200 ms. For the 
example shown in figure |^ we obtain an initiation time 
of ti = (215 ± 10) ms. This time corresponds to 30 tq, 
where tq — {nav)~^ is the classical collision time in a 
thermal gas (n: mean atomic density, a: scattering cross 
section, v: average thermal velocity). For the homoge- 
neous Bose gas the first stage of formation was predicted 
to be on the order of tq Q. Early numerical calcula- 
tions have confirmed that this should also be valid for the 
case of a harmonic confining potential, provided that all 
scattering processes are included [@,||. Most theoretical 
treatments require sufficient ergodicity, meaning that er- 
godic mixing occurs on the order of a few collision times. 
At temperatures below a few /iK this does not describe 
the experimental situation. Evaporation in the magnetic 
trap becomes one-dimensional due to the gravitational 
sag of the atom cloud. This leads to a significantly lower 
evaporation rate and the ergodic mixing times can 
be on the order of several ten collision times |10|. 
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The measured initiation times show a pronounced di- 
vergence at 77 ~ 4.6 and become much longer than the 
ergodic mixing time of some ten tq. Very slow cooling 
through the phase transition does not only affect the ini- 
tiation time but it also changes the characteristics of the 
growth curve. We detected a slow, approximately linear 
growth of a condensate prior to the initiation time, de- 
termined from equation (j^) . The open circles in figure ^ 
indicate when a condensate is observed for the first time 
and the gray stars show the time the condensate popula- 
tion exceeds 10^ as obtained from theoretical simulations 



of the growth experiments |14|. 

In figure ^ the two distinct stages of condensate for- 
mation can be seen for the case of slow cooling (77 = 4.6): 
first, there is a approximately linear growth at constant 
rate of Nq — (44 ± 8) x 10^ s~-^, followed by a very rapid 
growth with an initial rate 8 times larger. The slow lin- 
ear growth of the condensate indicates a less pronounced 
bosonic stimulation during this stage of the condensate 
formation. The numerical calculation does not reproduce 
the measured condensate growth data accurately in this 
regime. The calculation rather shows a growth rate of 
the same order of magnitude as the second stage of the 
observed growth, suggesting that the model works well 
if the system is sufficiently far below the transition tem- 
perature. For the thermal cloud the calculated temper- 
atures agree with the experimentally determined values. 
We deduce that a possible heating rate of the atoms in 
the trap is below our sensitivity of 5 nK/s and thus prob- 
ably negligible. The measured decrease in the total atom 
number for long time scales seems to be slightly larger 
than expected from the simulations. Whereas three body 
losses in the condensate are included in the calculations, 
we neglect three body recombination involving thermal 
atoms. The factor of 6 enhancement for these processes 
is negated by the much lower density of the thermal cloud 
cloud: for the given trap parameters and a temperature 
of 300 uK one obtains a lifetime of more than 40 s 
Similarly, the background limited lifetime of the magnetic 
trap is at least 40 s. 

We attribute the discrepancy between the theoretical 
model and the measured condensate numbers for slow 
cooling to approximations in the theoretical model which 
for example do not take into account coherence effects in 
the scattering process between the atoms. Possibly, the 
two-stage growth can be explained by the formation of 
quasi-condensates which have been predicted as an early 
stage in the formation process of a Bose-Einstein con- 
densate 1^] . Quasi-condensates exhibit strong phase fluc- 
tuations that die out as the system grows. These phase 
fluctuations are due to a population of low-lying excited 
states in the trap. Density fluctuations are suppressed 
due to their comparatively large excitation energy and 
therefore the density exhibits a Thomas-Fermi distribu- 
tion. The growth rate could then be significantly smaller 
than expected from the assumption of Bose-stimulation. 



Indeed, we find the growth rate to be almost one or- 
der of magnitude smaller than calculated from our model 
based on quantum kinetic theory. The time scale for both 
stages of the evolution is longer than the ergodic mixing 
time. The system is closer to equilibrium compared to 
the case of rapid cooling, i.e. small values of the evapora- 
tion parameter 77. However, the kink in the growth curves 
clearly shows that the condensate is not formed adiabat- 
ically. In this case the atom number should increase ac- 
cording to Nq = A^(l — {T/TcY) when the temperature 
is lowered smoothly. 

It was recently found that 3D Bose-Einstein conden- 
sates can exhibit strong phase fluctuations in elongated 
traps ijl^. The phase fluctuations are characterized by 
the parameter 5\ = {T /Tc)[N/Nof/^5l, where 5l is a 
factor that mainly depends on the trap geometry and 
atom number. Its value is about 0.15 for our exper- 
iment. The kink in the growth curves is observed at 
temperatures higher than 0.9 Tc, thus large phase fluctu- 
ations 5\ > \ can be expected in the regime Nq/N < 
4%. We flnd the transition from linear growth to re- 
laxational growth for slow cooling always in the range 
1.5% < No/N < 3.5%, which is compatible with the the- 
oretical expectation. When the phase fluctuations have 
died out, i.e. <C 1, a true Bose-Einstein condensate 
forms. This stage of the evolution then shows the same 
growth characteristics as observed for rapid cooling. So 
far calculation of the kinetics of the growth process have 
not taken into account the large aspect ratio of cigar 
shaped traps. Therefore axial phase fluctuations may be 
not included in these models. Evidence for axial phase 
fluctuations in a steady state of a very elongated conden- 
sate is reported in pO[ |. 

In conclusion, we have investigated the influence of the 
thermal atom cloud on the growth of a Bose-Einstein 
condensate. The measured data has been compared to a 
numerical model based on quantum kinetic theory. We 
find the minimum initiation time for condensate growth 
to be 30 times the classical collision time, in quantita- 
tive agreement with our calculations. This comparatively 
long timescale is presumably due to slow ergodic mixing 
of the trapped atomic cloud. Very slow cooling of the 
cloud enables us to observe a two stage process during 
the formation of a Bose-Einstein condensate, which is 
possibly beyond the approximations made in the model. 
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FIG. 1. Growth curve of a Bose-Einstein condensate for 
an evaporation parameter rj — 1.4. The lines are the re- 
sults of the numerical simulation of the growth for the start- 
ing conditions Ni = 4.4 x 10*^ and Ti = 610 nK (solid) and 
Ni = 4.2 X 10'' and Ti = 610 nK (dotted). Every data point is 
averaged over 3 identical repetitions of the experiment with 
statistical errors shown by the bars. 
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FIG. 2. Initiation times for the growth of Bose-Einstein 
condensates. The squares show the measured initiation times 
according to equation (P), the open circles show the times, 
when a condensate is detected for the first time, and the gray 
stars show the time the condensate population exceeds 10"* 
as obtained from theoretical simulations. The inset shows an 
enlargement of the results for small r/. 
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FIG. 3. Condensate growth curve for = 4.6. The kink 
in the growth curve is clearly seen. It occurs at a condensate 
fraction of No/N = 3.5%. The linear growth rate is inferred 
from a linear fit (straight line). The relaxational part of the 
growth curve is fitted by equation (^) (curved line). The 
initial rate in this regime is eight times larger than the linear 
growth rate. Every data point is averaged over 4 identical 
repetitions of the experiment with statistical errors shown by 
the bars. The results of the numerical calculation are shown 
for the same initial parameters as in fig. 0: A*; = 4.4 x 10^ and 
Ti = 610 nK (dashed) and Ni = 4.2 x 10*^ and Ti = 610 nK 
(dotted). 
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